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ABSTRACT 

An exotic dark component, named generalized Chaplygin gas (Cg) and pa- 
rameterized by an equation of state p = —A/p^^ where A and a are arbitrary 
constants, is one of the possible candidates for dark energy as well as for a uni- 
fied scenario of dark matter /energy. In this paper we investigate qualitative and 
quantitative aspects of the angular size - redshift test in cosmological models 
driven by such a dark component. We discuss the prospects for constraining 
the Cg equation of state from measurements of the angular size at low and high 
redshift radio sources and also from a joint analysis involving angular size and 
supernova data. A detailed discussion about the influence of the Cg on the mini- 
mal redshift at which the angular size of an extragalactic source takes its minimal 
value is also presented. 

Subject headings: Cosmology: Cosmological Parameters, Dark Matter - Equation 
of State 



1. Introduction 



The impressive convergence of recent observational facts along with some apparently 
successful theoretical predictions seem to indicate that the simple picture provided by the 
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standard cold dark matter model (SCDM) is insufficient to describe the present stage of 
our Universe. From these results, the most plausible picture for our world is a spatially 
flat scenario dominated basically by cold dark matter and an exotic component endowed 
with large negative pressure, usually named dark energy. Despite the good observational 
indications for the existence of these two components, their physical properties constitute a 
completely open question at present, which gives rise to the so-called dark matter and dark 
energy problems. 

Dark matter, whose the leading particle candidates are the axions and the neutralinos, 
was originally inferred from galactic rotation curves which show a general behavior that 
is significantly different from the one predicted by Newtonian mechanics. Dark energy or 
quintessence, whose the main candidates are a cosmological constant A and a relic scalar field 
0, has been inferred from a combination of astronomical observations which includes distance 
measurements of type la supernovae (SNe la) indicating that the Universe is speeding up 
not slowing down (Perlmutter et al. 1999; Riess et al. 1998; 2004), cosmic microwave 
background (CMB) data suggesting VLtotai — 1 (de Bernardis et al. 2000; Spergel et al. 2003) 
and clustering estimates providing Qm — 0.3 (Calberg et al. 1996; Dekel et al. 1997). While 
the combination of these two latter results implies the existence of a smooth component of 
energy that contributes with ~ 2/3 of the critical density, the SNe la results require this 
component to have a negative pressure thereby leading to a repulsive gravity (for reviews 
see Peebles & Ratra 2003, Padmanabhan 2003). 

By assuming the existence of these two dominant forms of energy in the Universe, 
one finds that the main distinction between them comes from their gravitational effects. 
Cold dark matter agglomerates at small scales whereas dark energy seems to be a smooth 
component. In certain sense, such properties are directly linked to the equation of state 
of both components (for a short review see Lima 2004). On the other hand, the idea of 
a unified description for the CDM and dark energy scenarios has received much attention 
(Matos & Urciia-Lopez 2000; Davidson, Karasik & Lederer 2001; Kasuya 2001; Watterich 
2002; Padmanabhan & Choudhury 2002). For example, Wcttcrich (2002) suggested that dark 
matter might consist of quintessence lumps while Kasuya (2001) showed that spintessence- 
like scenarios are generally unstable to formation of Q balls which behave as pressureless 
matter. More recently, Padmanabhan and Choudhury (2002) investigated such a possibility 
via a string theory motivated tachyonic field. 

Another interesting attempt of unification was suggested by Kamenshchik et al. (2001) 
and developed by Bilic et al. (2002) and Bento et al. (2002). It refers to an exotic fiuid, the 
so-called Chaplygin gas (Cg), whose equation of state is given by 



PCg = -A/p'^g, 



(1) 
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with a — 1 and A a positive constant related to the present day Chaplygin adiabatic sound 
speed, Vg = aA/pcg^o, where pcg,o is the current Cg density. In actual fact, the above 
equation for a ^ 1 constitutes a generalization of the original Chaplygin gas equation of 
state proposed by Bento et al. (2002). The idea of a dark- matter-energy unification from an 
equation of state like Eq. (1) comes from the fact that the Cg or the generalized Chaplygin 
gas (from now on we use Cg to denote the Chaplygin gas as well as the generalized Chaplygin 
gas) can naturally interpolate between nonrelativistic matter and negative-pressure dark 
energy regimes. It can be easily seen by inserting the Eq. (1) into the energy conservation 
law (li^T.r = 0). One finds, 

1 

l+a 

(2) 

where the subscript o denotes present day quantities, R{t) is the cosmological scale factor, 
and As = A/ p^^'^ is a quantity related to the present sound speed for the Chaplygin gas 
{Vg = a As). As can be seen from the above equations, the Chaplygin gas interpolates 
between non-relativistic matter [pcg{R — >^ 0) oc R~^] and negative-pressure dark component 
regimes [pcg{R oo) oc const.]. 

Motivated by these potentialities, there has been growing interest in exploring theoret- 
ical (Bordemann & Hoppe, 1993; Hoppe 1993; Jackiw 2000; Gonzalez-Diaz 2003a; 2003b; 
Krcmcr 2003; Khalatnikov 2003; Balakin et al. 2003; Bilic et al. 2003) and observational 
consequences of the Chaplygin gas, not only as a possibility of unification for dark mat- 
ter/energy but also as a new candidate for dark energy only. The viability of such scenarios 
has been tested by a number of cosmological tests, including SNe la data (Fabris, Goncalves 
& de Souza, 2002; Cohstete Jr. et al. 2003; Avehno et al. 2003; Makler, de Oliveira & 
Waga 2003; Bertolami et al. 2004), lensing statistics (Dev, Alcaniz & Jain 2003; Silva & 
Bertolami 2003; Dev, Alcaniz & Jain 2004), CMB measurements (Bento, Bertolami & Sen 
2003a; 2003b; 2003c; Carturan & Finelh 2002; Amendola et al. 2003), age-redshift test 
(Alcaniz, Jain & Dev 2002), measurements of X-ray luminosity of galaxy clusters (Cunha, 
Lima & Alcaniz 2003), future lensing and SNe la experiments (Avelino et al. 2003; Silva 
& Bertolami 2003; Sahni et al. 2003), statefinder parameters (Sahni et al. 2003), as well 
as by observations of large scale structure (Multamaki, Manera & Gaztanaga 2003; Bihc 
et al. 2003; Bega et al. 2003; Bean and Dore 2003; Avehno et al. 2004). In particular, 
the latter reference has shown that in the context of unified dark matter/energy models the 
onset of the non-linear regime on small cosmological scales may lead to the breakdown of 
the background solution, even on large cosmological scales (if a is not equal to zero). 

Although carefully investigated in many of its theoretical and observational aspects, 
the infiuence of a Cg component in some kinematic tests such as the angular size-redshift 



PCg - PCga 



3(l+a) 
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relation {9 — z) still remains to be studied. This is the goal of the present paper. In 
the next sections we investigate qualitative and quantitative aspects of the angular size- 
redshift test in cosmological models driven by this dark matter/energy component. We first 
investigate the influence of the Cg on the minimal redshift at which the angular size of an 
extragalactic source takes its minimal value. Afterwards, we consider the 9{z) data recently 
updated and extended by Gurvits et al. (1999) to constrain the equation of state of the 
Cg component as well as a combination between these 9{z) observations and the latest SNe 
la data, as provided by Riess et al. (2004). To do so, we follow Alcaniz et al. (2002) 
and consider two different cases, namely, a flat scenario in which the generalized Chaplygin 
gas together with the observed baryonic content are responsible by the dynamics of the 
present-day Universe [unifying dark matter-energy] (UDME) and a flat scenario driven by 
non-relativistic matter plus the generalized Chaplygin gas (CgCDM). For UDME scenarios 
we adopt in our computations Vt^, = 0.04, in accordance with the estimates of the baryon 
density at nucleosynthesis (Buries, Nollett & Turner 2001) and the latest measurements of 
the Hubble parameter (Freedman et al. 2002). For CgCDM models we assume Qm = 0.3, 
as suggested by dynamical estimates on scales up to about 2h~^ Mpc (Calberg et al. 1996; 
Dekel et al. 1997). For the sake of completeness an additional analysis for the original Cg 
model {a — 1) is also included. 

This paper is organized as follows. In Sec. II we present the basic equations and 
distance formulas necessary for our analysis. Following the method developed in Lima & 
Alcaniz (2000a; 2000b) the influence of a Cg-like component on the minimal redshift Zm 
is investigated. In Sec. Ill we analyze the constraints on the equation of state of the Cg 
component from measurements of the angular size of compact radio sources and SNe la data 
and compare them with other independent limits. We end the paper by summarizing the 
main results in the conclusion Section. 



2. 



Distance Formulas and the minimal redshift 



2.1. 



Distance Formulas 



Let us now consider the flat FRW line element (c = 1) 



(3) 



where Xi ^"^^ ^-^e dimensionless comoving coordinates. In this background, the angular 
size-redshift relation for a rod of intrinsic length £ is easily obtained by integrating the spatial 
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part of the above expression for x ^-nd fixed, i.e., 



= -. (4) 



The angular diameter distance (Ia is given by 

^ (1 + ^) {l + z)J^,^,^-.E{n„A,a,x) ' ^ > 



where x — = (1 + z)"^ is a convenient integration variable. For the kind of models here 
considered, the dimensionless function E{D,j, Ag, a, x) takes the following form 

^ 1 >t 1/2 

o.) A 

llj) ^3(a+l) 




(6) 



where Qj stands for the baryonic matter density parameter (j = b) in UDME scenarios and 
the baryonic + dark matter density parameter (j — m) in CgCDM models. Note that from 
the above equation UDME models reduce to the ACDM case for a = whereas CgCDM 
models get the same hmit, regardless of the value of a, when Ag — 1. In both cases, the 
standard Einstein-de Sitter behavior is fully recovered for Ag — Q (Avelino et al. 2003b; 
Fabris et al. 2004). 



2.2. Minimal redshift 

As widely known, the existence of a minimal redshift on the angular size-redshift 
relation may qualitatively be understood in terms of an expanding space: the light observed 
today from a source at high z was emitted when the object was closer (for a pedagogical 
review on this topic see Janis 1986). The relevant aspect here is to demonstrate how this 
effect may be quantified in terms of the Chaplygin parameters a and As. To analyze the 
scnsivity of the critical redshift to this dark component, we adopt here the approach originally 
presented by (Lima & Alcaniz 2000a; 2000b). The numerical results of this method have 
been confirmed by Lewis & Ibata (2002) through a Monte Carlo analysis. 

The redshift z„i at which the angular size takes its minimal value is the one cancelling 
out the derivative of 9 with respect to z. Hence, from Eqs. (4)- (6) we have the following 
condition 

Xm (1 ~l~ ^m)X. \z=Zm.i C^) 

where a prime denotes differentiation with respect to z and by definition Xm = xi^m)- Note 
that Eq. (5) can readily be differentiated, yielding 

(1 + Zm)x'\z=zm = -B-FT^i^j^ ^s, a, Zm), (8) 
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where 



^ = { 



(1 + Zm) 



-3 + 



(1 - As) 



a + l 




(9) 



(l + zJ-3("+l)J 



Finally, by combining equations (7)- (9), we find 




(^1 Zm) 



dx 



(10) 



In Fig. 1 we show the results of the above expression. Panels (a) and (c) show the 
minimal redshift function of the parameter Ag for selected values of a in the context 

of CgCDM {VLm = 0.3) and UDME (fi^ = 0.04) models, respectively, whereas Panels (b) and 
(d) display the a — Zm plane for some values of Ag also for CgCDM and UDME scenarios 
(note that the scale of Zm in each panel is different). The minimal redshift is a much more 
sensitive function to the parameter ^1^ than to the index a. As can be seen in the panels [and 
also expected from Eq. (5)], regardless of the value of a, models with Ag — reduce to the 
Einstein-de Sitter case so that the standard result Zm — 1-25 is fully recovered. As physically 
expected (see, e.g., Lima & Alcaniz 2000a; 2000b), the smaller the contribution of the 
material component (baryonic and/or dark) the higher the minimal redshift z^- From this 
qualitative argument UDME models would be in a better agreement with the observational 
data than CgCDM scenarios since the current data for milliarcsecond radio sources do not 
show clear evidence for a minimal angular size (6'minimai) close to z — 1.25 (Gurvits et al. 
1999). For the best fit CgCDM model obtained from an analysis involving galaxy clusters 
X-ray and supernova data, i.e., Ag — 0.98 and a — 0.93 we find Zm — 1-6, a value that is 
very similar to the one predicted by the current concordance model, namely, a flat ACDM 
scenario with fl^ = 0.3 and also by the standard FRW model with flm = 0.5 (see Table 
I of Lima & Alcaniz 2000a). It is worth mentioning that if high-redshift sources present 
cosmological evolution, such an effect would move the position of the ^minimal to extremely 
high redshift . 



In this section we study the constraints from angular size measurements of high- 2; radio 
sources on the free parameters of the Cg model. In order to constrain the parameters Ag 
and a we use the angular size data for milliarcsecond radio sources recently compiled by 



3. Constraints from angular size measurements 
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Gurvits et al. (1999). This data set is composed by 145 sources at low and high redshifts 
(0.011 < z < 4.72) distributed into 12 bins with 12-13 sources per bin. Since the main 
difference between the analysis performed in this Section and the previous ones that have used 
these angular size data to constrain cosmological parameters is the background cosmology, 
we refer the reader to previous works for a more complete analysis and detailed formulas 
(see, for instance, Jackson & Dodgson 1997; 2003; Vishwakarma 2001; Lima & Alcaniz 2002; 
Zhu & Fujimoto 2002; Jain et al. 2003; Chen & Ratra 2003). 

Following a procedure similar to that described in the quoted reference, we determine 
the cosmological parameters As and a through a minimization for a range of As and a 
spanning the interval [0, 1] in steps of 0.02, i.e.. 



where 9{zi, Vtj, As, a) is given by Eqs. (4)-(6) and Oobsi is the observed values of the angular 
size with errors CTj of the ith bin in the sample. 68% and 95% confidence regions are defined 
by the conventional two-parameters levels 2.30 and 6.17, respectively. 

Figures 2a and 2b show the binned data of the median angular size plotted as a function 
of redshift for selected values of Ag and a for UDME and CgCDM models, respectively. For 
comparison the current "concordance" scenario, i.e, the fiat ACDM model with 0,^ — 0.3 is 
also shown (thick hne). In Figs. 3a and 3b, we show the 68% and 95% confidence contours 
computed by fixing the characteristic length D at its best fit value from a minimization of 
Eq. (11) relative to the parameters Ag, a and D. From this analysis we obtain the following 
best fit values: As = 0.96, a = 0.84 and D = 2A.lh-^ pc (xLn = 4.44) for UDME scenarios 
and As = 1.0 and D = 29.58/i~^ pc (Xmm = 4.57) for CgCDM models. Both cases represent 
accelerating universes with deceleration parameter and the total age of the Universe given by 
Qo = -0.88 and to ~ llh'^ Gyr (UDME) and Qo = -0.55 and to ~ QAh'^ Gyr (CgCDM). 
Note that, as happens in the calculation of the minimal redshift Zm as well as in other 
cosmological tests (see, e.g., Cunha, Lima & Alcaniz 2004), the current angular size data 
constrain more strongly the parameter Ag than the index a. From Fig. 3 it is perceptible 
that while the parameter As is constrained to be > 0.68 (UDME) and > 0.54 (CgCDM) at 
2(7 the entire interval of a is allowed. 



By combining the angular and luminosity distances, interesting constraints on the Cg 
parameters are obtained. To perform such analysis, we follow the conventional magnitude- 




(11) 



3.1. 9{z) + SNe la 
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redshift test (see, e.g, Goliath et al. 2001; Dicus & Repko 2003; Alcaniz & Pires 2004) and 

use the latest SNe la data set that corresponds to the gold sample (157 events including 9 
SNe at 2; > 1 ) of Riess et al. (2004). In this analysis, both the characteristic angular size D 
and the Hubble parameter are considered "nuisance" parameters so that we marginalize 
over them. In the case of the d{z) test the marginalization over the characteristic length can 
be easily done by defining a modified statistics as: 



where 



and 




For a similar procedure concerning the magnitude-redshift test we refer the reader to (Gohath 
et al. 2001; Bertolami & Silva 2003). 

Figures (4) and (5) show the results of our analysis. In Panels 4a and 4b we display 
contours of the combined likelihood analysis for the parametric space a — Agva. the context 
of CgCDM and UDME scenarios, respectively. For UDME models we see that the available 
parameter space is considerably modified when compared with Fig. 3b, with the best-fit 
value for Ag provided by the d{z) analysis, i.e., Ag = 0.96 being off by ~ 3a relative to 
the joint analysis. This analysis also yields 0.65 < Ag < 0.90 (95% c.l.) for UDME models 
and As > 0.85 (95% c.l.) for CgCDM scenarios. These particular limits on As are in good 
agreement with the ones obtained from quasar lensing statistics {As > 0.72, Dev, Jain & 
Alcaniz 2002), the old SNe la data (A, = 0.87t||;}^, Fabris et al. 2002; Avelino et al. 2002) 
as well as from the location of the acoustic peaks of CMB as given by BOOMERANG 
and Arqueops (0.57 < Ag < 0.91 for a < 1, Bento et al. 2003). However, they are only 
marginally compatible with the tight constraint obtained from the expected number of lensed 
radio source for the Cosmic All-Sky Survey (CLASS) statistical data (Dev, Jain & Alcaniz 
2004). Other interesting limits from this analysis are also obtained on the original version of 
Cg, i.e., by fixing a = 1. In this case, the plane — Ag (Fig. 5a) is reasonably constrained 
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with the best fit values located at Ag — 0.81 and = 0.0 with xLin/^ — This 
particular value of the matter density parameter agrees with the one found by Fabris et 
al. (2002) by using the old sample of SNe la from High-z Supernova Project. In this Ref., 
the authors interpret such a result as a possible backup to the idea of dark matter-energy 
unification (UDME models). However, in the light of recent CMB data, unification from the 
original Cg (a = 1) seems to be quite unrealistic since the location of the acoustic peaks as 
given by WMAP and BOOMERANG is in conflict with the predictions of this particular 
scenario (see, e.g., Carturan & FineUi 2003). Figure 5b shows the same analysis of Panel 
a by assuming the Gaussian prior on the matter density parameter, i.e., — 0.27 ± 0.04 
(Spergel et al. 2003). As can be seen, the parameter space is tightly constrained with the 
best-fit scenario located at Vtm = 0.26, Ag = 0.97 and Xmin/^ — l-l- Table I we present 
the best fit values for the Cg parameters obtained from the two main analyses performed in 
this paper. 



4. Conclusion 



Based on a large body of observational evidence, a consensus is beginning to emerge 
that we live in a fiat, accelerated universe composed of ~ 1/3 of matter (barionic -|- dark) 
and ~ 2/3 of a negative-pressure dark component. However, since the nature of these dark 
components (matter and energy) is not well understood, an important task nowadays in 
cosmology is to investigate the existing possibilities in the light of the current observational 
data. In this paper we have focused our attention on some observational aspects of cosmolo- 
gies driven by an exotic dark energy component named generalized Chaplygin gas. These 
models also constitute an interesting possibility of unification for dark matter/energy (where 
these two dark components are seen as different manifestations of a single fluid). Initially, 
We have investigated the influence of the Cg on the minimal redshift (zm) at which the 
angular sizes of extragalactic sources takes their minimal values. From this analysis it was 
showed that the location of the minimal redshift is a much more sensitive function to the 
parameter Ag than to the index a. By using a large sample of milliarcsecond radio sources 
recently updated and extended by Gurvits et al. (1999) along with the latest SNe la data, 
as given by Riess et al. (2004) we obtained, as the best flt for these data, Ag — 0.84 and 
a = 1.0 (UDME) and Ag = 0.99 and a = 1.0 (CgCDM). Such values are fully disagreement 
with the CMB analysis by Amendola et al. (2004) which showed that Cg scenarios with 
a > 0.2 are ruled out by the current WMAP data. Finally, it should be remarked that the 
background results presented here may be somewhat modifled if one takes into account the 
onset of the non- linear regime, as recently discussed by Avelino et al. (2004). 
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Fig. 1. — Critical redshift function of the Cg parameters Ag and a for CgCDM 

(Panels a and b) and UDME (Panels c and d) scenarios. Note that, regardless of the value 
of a, all models with = recover the standard Einstein-de Sitter result, i.e., Zm — 1-25. 
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Fig. 2. — Angular size versus redshift for some combinations of the parameters As and a 
for CgCDM and UDME models. The characteristic length D has been fixed at 29.58/i~^ pc 
(Panel a) and at 24.10/i^^ pc (Panel b), which corresponds to the best fit values obtained 
from a minimization of Eq. (11) relative to the parameters Ag, a and D. 
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Fig. 3. — Confidence regions in the a — Ag plane according to the updated sample of angular 
size data of Gurvits et al. (1999) for CgCDM (Panel a) and UDME (Panel b) models. From 
this figure we see that while the entire interval of a is allowed the parameter As is constrained 
to be > 0.54 (CgCDM) and > 0.68 (UDME) at 95% confidence level. 
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Fig. 4. — a) 68% and 95% confidence levels in the plane a — As for CgCDM models corre- 
sponding to the joint 6{z) + SNe analysis as described in the text. As in Fig. 3, the entire 
interval of a is allowed. The parameter Ag is now constrained to be > 0.85 at 95% c.l. b) 
The same as in Panel a for UDME scenarios. For this case we find 0.65 < Ag < 0.90 (95% 
c.l.) 
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Fig. 5. — a) The Q,n ~ plane for the joint 9{z) + SNe analysis having a = 1 (original 
version of Cg). Here the best fit values are located at Ag = 0.81 and = 0.0. b) The same 
as in Panel a by assuming a Gaussian prior on the matter density parameter flj^ = 0.27±0.04 
(Spergel et al. 2003). The best-fit scenario is located at Q.m ~ 0.26 and Ag — 0.97 with 
• /i/ ~ 1.1. 
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Table 1: Best fit values for Cg parameters 
Test UDME CgCDM 



9{z) ^ = 0.96 a = 0.84 ^ = 1.0 a = 0.96 

SNe A = 0.84 a = 1.0 = 0.98 a = 0.96 

9{z) + SNe A, = 0.84 a = 1.0 A^ = 0.99 a = 1.0 



